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We discuss the possible role of isocurvature perturbations for the quantum decoherence of the 
curvature perturbation during inflation. We point out that if the inflaton trajectory in field space is 
curved, the adiabatic mode is generically coupled to the isocurvature mode and thus tracing out the 
latter can cause the curvature perturbation to decohere. We explicitly investigate this suggestion 
in a model of inflation with two decoupled massive fields and find that decoherence is effective for a 
wide range of mass ratios. In particular, we calculate the entanglement entropy for this model and 
show that it grows after horizon exit, providing a quantitative measure for decoherence. 
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Inflation has become the leading paradigm of early universe cosmology not least because of its ability to imprint 
scale invariant inhomogeneities on superhorizon scales via a causal mechanism [l[ . These inhomogeneities are thought 
to provide the seeds which later become the temperature anisotropics in the Cosmic Microwave Background and the 
Large Scale Structure in the Universe. In inflation, the origin of these initial perturbations lies in quantum fluctuations 
' of matter fields, amplified during a period of quasi-exponential expansion, which source adiabatic and isocurvature 
cosmological perturbations. The adiabatic perturbation is a variation in total energy density or equivalently in the 
gravitational potential, while isocurvature perturbations are variations in the entropy density of various matter com- 
ponents. All inflationary models share this basic mechanism, with different models characterized by small deviations 
from exact scale invariance. Along with other calculable properties, these deviations are particularly appealing as 
they allow for inflationary models to be tested against increasingly precise observations 0- 
. Inflationary calculations are based on the quantum mechanics of scalar fields in expanding spacetimes [||, where 
t— I ' the relevant observable is the amplitude of the field's Fourier modes. Although treated as a quantum mechanical 
variable, this amplitude is interpreted as a stochastic random variable described by a gaussian distribution, with the 
variance given by the power-spectrum. This interpretation, used in CMB analyses and simulations of Large Scale 
Structure, proves to be very accurate for calculational purposes; quantum correlation functions, exactly calculated, 
fS are practically the same as the corresponding correlation functions obtained by using the relevant classical probability 
distribution (4J. 

The consistency of this stochastic interpretation requires a density matrix which is diagonal in the amplitude basis. 
However, the density matrix of inflationary perturbations is not automatically diagonal in this basis (see eg. [5[). 
One cannot assign a specific amplitude to the Fourier modes which must thus be considered to exist in coherent 
superpositions of different amplitudes. Prior to a stochastic interpretation, an as-of-yet undisclosed decoherence 
. — i , mechanism is required to diagonalize the density matrix. 

The phenomenon of decoherence is ubiquitous in quantum physics and originates from couplings of the system 
of interest with degrees of freedom belonging to some unobservable environment d 0|- Various arguments and 
calculations suggesting that a form of such environmental decoherence can indeed occur for inflationary perturbations 
have been put forward in @, ©, [H [U 03, 0, El El El • In [1, El El El El El El > the environment is taken to 
consist of the short wavelength modes which are coupled to the long wavelength modes via non-linear couplings. In 
this paper, we point out that the decoherence of cosmological perturbations can occur in a precisely defined dynamical 
setting if the role of the environment is played by unobservable isocurvature perturbations generated during multi-field 
inflation. The coupling in this case are bilinear and no non-linearities are invoked. 
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II. QUANTUM MECHANICS OF INFLATIONARY PERTURBATIONS 



To study decoherence in inflation it is convenient to use the Schrodinger picture which we briefly review in this 
section. The action for a free massive scalar field is 



m 2 (j) 2 ) 



(1) 
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where we take the spacetime metric to be 

ds 2 = g^dx^dx" = -dt 2 + a{t) 2 5 ll dx l dx j . (2) 
Then, the canonical momentum conjugate to the field 0, 

i)C 



allows us to write the Hamiltonian H = J d 3 x (ivdt4> — C) 



= a d d t 4>, (3) 



H = Jd 3 x^ (^+a(d t 



+ a 3 m 2 c/) 2 . (4) 



Quantization dictates the replacement 7r(x) — > — ih which realizes the commutation relation [</>(x), 7r(y)] = 
ih5(x. — y). The wave functional \&[0] obeys the Schrodinger equation 

d 

ih—y = H^. (5) 

dt ' 

Since we will be considering perturbations in their ground states, we take the wavefunctional to have a gaussian form 



= TVexp 



^2;cPi/0(x),4(x,y,i)0(y) 



(G) 



with M the normalization factor. Due to the homogeneity of the vacuum state, the correlator in ^ satisfies, 
A(x, y, t) = A(y,x, t). For this state the Schrodinger equation gives 

h 2 r 

ihdtlnM = — / d 3 xA(x,K,t) (7) 

h 2 r 

ihd t A(y, z,t) = — d 3 xA(x, y, t)A(x, z,t)+a (d 2 - a 2 m 2 ) <5 (3) (y - z) . (8) 

Since the theory is noninteracting, it is convenient to solve (|8]) in momentum space, where it becomes local. Indeed, 
upon writing 

d 3 k 



(x) = J j^ste**, (9) 

/d 3 k 
— A(k,ty k -^-y\ (io) 



with = 0k* and A(— k, t) = A(k, t), we get 

h 2 



ihd t A(k,t) = -^A 2 {\i,t) - a{k 2 + a 2 m 2 ) . (11) 



The Heisenberg picture mode functions ipk(r]) of a scalar field with the action |T]) satisfy 

d 2 (aMv)) + (k 2 + ^r^-) WW) = o , (12) 

where r\ is the conformal time, defined via adr\ = dt, and v 2 — j — jjj. The Bunch-Davies vacuum solution of (fT2l is 
of the form 

aMv) = ] f^e^ + ^Hi 1 \~k V ) , (13) 
and the correlator A(k) (fTTj) is related to the mode functions (fl"3"| via [|[ 

A(k ^ ) = 27# k w( 1 ^ a ^ l ^ )|2 )- (14) 



3 



At early times, in the regime k 2 n 2 3> 1 when the mode is "deep inside the horizon" and spacetime curvature is not 
important, the solution (fl"3"|) or equivalently ^ and p4[) reduces to that of Minkowski vacuum. At late times, when 
k 2 T] 2 <§; 1, the mode has "exited the horizon" and amplification ("particle production") occurs. The wave function 
for each mode can be written as 1 

*(<?!) k ,77) cxexp ^-i</. k v4(k,77)0*^ , (15) 

where the normalization factor can be easily obtained from ([7]) and by requiring J dtfidcf)* | "J | 2 = 1. From now on 
we shall be ignoring such factors. Note that in the Schrodinger picture the time evolution of the wave function is 
determined by the correlator A(k, n); k is time independent and should not be confused with the Heisenberg picture 
mode Vk(f?) (fT^I - 

The quantum mechanics of cosmological perturbations has been studied in (3), [13] • The system evolves into a kind 
of state known as a squeezed state, which exhibits a high degree of WKB classicality 

d q S(q)V, (16) 

where q and p are the configuration variable and its conjugate momentum, while S is the exponent of the wave function, 
practically the classical action. Thus, the amplitude of each mode and its conjugate momentum are related to a very 
high degree of accuracy via the corresponding classical relation. In cosmology, this state of affairs is interpreted as 
equivalent to a statistical mixture of mode amplitude eigenstates, where the probability that any of these states has 
been realized in our universe is given by ^(t^, ij)\ 2 . This would correspond to a density matrix which is diagonal in 
the field amplitude basis 

p k = ^|*(^ k ,r 7 )| 2 |^ k )^k|. (17) 

Such a stochastic mixture is a very good approximation when calculating correlators such as (0 k ) 4], which are then 
related to the stochastic properties of cosmological perturbations. However, the wave function (fTS"]) does not directly 
lead to the interpretation suggested by (fT"7|) . Explicitly calculating the density matrix from (fT"5|) , one finds 

(0k|p k |0k) = p(0k,0k) = *(0 k ,77)**(0k,77) oc exp (- u k ^[A]ui- iA k K[A]A£ - i3L4Mu k A k ]), 

where we have defined u = (cf> + 4>)/2 and A = <p — <f> . The above expression does not vanish for A ^ 0, making 
the off-diagonal terms of the density matrix non-zero. Thus, strictly speaking, the description (|17[) is not valid and 
one cannot say that any particular mode amplitude has been realized with a certain probability. The modes of 
cosmological perturbations exist in coherent superpositions with different mode amplitudes. In order to diagonalizc 
the density matrix, some process of decoherence must take place. As we show below, the existence of fields other 
than the inflaton during inflation can decohere the density matrix of cosmological perturbations even if no direct 
coupling between the different fields is assumed. The relevant coupling is ensured by the existence of scalar metric 
perturbations. 

III. ISOCURVATURE PERTURBATIONS AND DECOHERENCE 

A. Adiabatic and isocurvature perturbations 

The existence of isocurvature modes is ubiquitous when more than one scalar fields are relevant during inflation. 
We will therefore consider a two-field model of inflation with a potential V((p,x)- For cosmological perturbations for 
the metric and the scalar fields 

0(x,t)=^(i)+^(x,t) ) X(x,t)=x(*)+*X(M), (18) 
ds 2 = -(1 + 2$(x, t))dt 2 + a 2 (t)(l - 2$(x, t))dx 2 , 



The reality of the field </>(x) implies that the modes </>_k and <frk are related as 0_k = k , which necessitates that wave functions be 
represented in terms of two-mode states Il7|l . 



4 



the Lagrangian describing the dynamics of gauge invariant perturbations can be written as fl8j 

L = Jd 3 x ifl^qft,q - ~q (-V 2 + (aH) 2 n) q , (19) 

where 

q = a{5(p + jj$)e v + a ( 6 X + ^ $ ) § x = Q^v + q x e x , (20) 
with e v and e x the unit basis vectors in the ip - x held coordinates. The matrix f2 is given by 

n = a J 2 e a ® e 6 - (2 - e)l - 2e ((3 + e) e x <g) ei + 2?? ll ei <g) e x + 77 1 - (ei <g) e 2 + e 2 <8> ei)J , (21) 

where e = —H/H 2 and 

11 o 9V ■ ei ± _ <9V • e 2 

77 11 = -3 : , 77 = : , (22) 

H\</>\ H\cj>\ 

arc two slow roll parameters related to the acceleration parallel and perpendicular to the field velocity [l|[. The 
subscripts {a, b} take the values {p, x} an d we have used the unit vectors ei and e 2 to be defined shortly. Note 
that even if scalar fields are decoupled and the potential term is diagonal, the q a fields are still coupled due to the 
inclusion of the metric perturbation <f>. Thus, perturbations are always coupled gravitationally. As we discuss below 
this coupling can suffice to decohere the perturbations if one degree of freedom is traced out. 

In order to proceed, it is convenient to define two new basis vectors, tangential and normal to the field trajectory 

ei = cos e v + sing e x = —===e v + X e x (23) 

vr+r Vr+r 

e 2 = — sin0e v + cos0e x , (24) 
and express the Lagrangian in terms of components in this basis - we have 

L = Jd 3 x\ (9,q + Zq) T (a„q + Zq) - iq T (- V 2 + (aH) 2 Q) q , (25) 

where (to, n = 1, 2) 

qi ) , q m = e m • q, (26) 
Q2J 

^mn — (^7) 

is the matrix l~i expressed in the {ei , e 2 } basis and 

Z = d v 6 (? ~S] . (28) 



1 / 

Then, using the conjugate momentum 7r(x) = g^^x)) the Hamiltonian can be written as 

H = Jd 3 x Vtt - Vzq - ±q T Z T Tr + iq T (-V 2 + (aH) 2 Q) q . (29) 

The quantum Hamiltonian is obtained via the replacement 7r(x) — > ~~i^ sqfx) anc ^ ^he wav sfunction satisfies the 
Schrodinger equation 

iW> = , (30) 

which, along with the vacuum ansatz 



* = A^cxp f-1 Jd 3 xd 3 y q T (x)B(x- y)q(y)j , 



(31) 
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leads to 

»fifyB(x-y) = h 2 (^y"d 3 zB(x-z)B(z-y)j +ift[B(x-y), Z] + (V 2 - (aff) 2 Sl) <5^ (x - y) , (32) 

ihd v lnN = | / ^B(O). (33) 

The solution to these equations provides the full quantum description of cosmological perturbations in the two-field 
system. 

B. The reduced density matrix 

The process of decoherence relies on the existence of degrees of freedom which couple to - get entangled with - the 
system of interest but which are not directly accessible. It is assumed that they cannot be measured along with the 
system and thus provide an uncontrolable environment. Tracing them out leads to the diagonalization of the density 
matrix for the system of interest and thus explains the appearance of a classical world where different alternatives - 
like Schrodinger's cat being dead and alive - do not interfere and are perceived as distinct [1, 0|- In Cosmology, the 
role of the cat is played by the long wavelength fluctuations which lead to structure formation. 

In what follows we suggest that decoherence can occur for cosmological perturbations if the isocurvature component 
<?2 is considered to be the unobservable environment and is traced out. In cosmological considerations the two degrees 
of freedom in q describing perturbations are not directly observable. A more tangible remnant of inflation is the 
curvature perturbation which persists in later epochs of cosmic evolution and is the reason for structure formation. 
More precisely, it is given by a scalar quantity TZ which parameterizes the spatial curvature of comoving time slices 
(see e.g. [13|). During inflation it is related to qi by 

U=^f=^, (34) 

Vv 2 + x 2 a 

and for single field inflation it is constant on long weavelengths. The presence of an isocurvature perturbation makes 
TZ evolve in time [H, [H[ 



;K = 4 ^ . (35) 



d ^ = 2 dV ■ e 2 £2 

dN ip 2 + x 2 a 

The hat acts as a reminder that we are dealing with quantum operators. 

Any process which couples to TZ during inflation is effectively measuring qi , the adiabatic component of the per- 
turbations. A measurement of g 2 via gravitational interactions alone would entail a coupling to -j^TZ during inflation 

as can be seen from (|35[) . Thus, the gravitational perturbation TZ, which does not have its own independent dy- 
namics, acts as an apparatus which "measures" - gets entangled with - q\. Let us now assume that after inflation 
both <j> and x decay into relativistic particles, generating no entropy perturbation and making TZ the only remnant 
from the inflationary epoch. 2 In this case, information about -rkfTZ during inflation is lost and all post-inflationary 

perturbations can be obtained from TZ at the end of inflation, which is conserved on super-Hubble scales after q 2 
decays. Once perturbations reenter the horizon at late times, -j^TZ becomes different from zero again. This late time 

behaviour reflects late time dynamics of matter fields on small scales and has nothing to do with the inflationary -j^TZ. 
Consequently, as information about qi is lost, we can consider it to be an inaccessible environment which cannot be 
directly observed and can therefore be traced out of the density matrix. 

Under the above assumption, the density matrix relevant for the adiabatic perturbation alone is the one obtained 
by tracing out q 2 

= / 'dQ5jdt?2*(?i>92)**(9i)52) , (36) 

where the wavefunction \I/ is given by (|3ip . Expressing the result in terms of the variables u = (qi + <Zi)/2 and 
A = q\ — q\ , we find 



p(u, A) = exp 



(37) 



2 This reasoning is not valid if a gravitational potential is generated after inflation via a curvaton type mechanism l'2 ll . |22H . 
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with 

c= ( *{mi]-?g$) i^-^r 1 ) 

"H^lBn]-S»d) * (*[ Bll ] + %?) 

As we remarked above, decoherence would imply that the A components are suppressed when |A| ^ 0. This will 
happen if 

Cn/C 22 «l, (39) 

since in this case the relative importance of the A over the u terms will diminish, making quantum correlations 
unimportant. Thus, we expect that decoherence will become effective if 

$t[B u ] > 2 v / &[B u ]ft[B 22 ] . (40) 

In inflationary models with two scalar fields, f2i 2 in (f3"2| is in general non zero, Vii 2 ^ and thus Bi 2 ^ (see also 
(|43p and (|46[) ) . In particular, for inflationary perturbations and on superhorizon scales, the imaginary part Q ! [-Bi 2 ] 
strongly dominates over the real parts 3?[Bn] and 5ft[-B 22 ] and thus, when isocurvature modes are inaccessible, we 
expect that decoherence will occur generically. 3 We check this assertion explicitly in the next section for a simple two 
field model with two massive free fields. We find that this expectation is indeed correct for a relatively wide range of 
parameter values of the model. 




IV. DECOHERENCE IN THE MODEL V(tp,x) = fmV + §mV 

In this section we illustrate the role of isocurvature perturbations for decoherence in an inflationary model with 
V((f, x) = \rn 2 ip 2 + \y?x 2 i a model which has been studied extensively in [23[ (see also [l8j]). We first present exact 
numerical results for two cases: a) /i = 1.05m and b) fj, = IOto. We then derive an analytic result for the first case. 
We note that the first case is practically the same as single-field inflation as the trajectory in field space is practically 
straight. In the second case, the heavier field drops out of slow roll and start oscillating approximately 40 e-folds 
before inflation ends. However, as cosmologically observable scales exit the horizon while both fields are still in slow 
roll, the spectral index does not deviate significantly from unity and thus the model is observationally viable. 

The wavefunction obtained from (|32|) explicitly reads: 



ihdr,B n = h 2 (Bfi + B? 2 ) + 2ihd, 1 6B 12 - (k 2 + (aH) 2 Q n ) , (41) 
ihd v B 22 = h 2 (B 2 2 + B? a ) - 2ihd v 8B 12 - (k 2 + (aH) 2 Q 22 ) , (42) 
ihd v B 12 = fi 2 B 12 (B 11 + B 22 )+^(B 22 -B 11 )a i; 0-(ai/) 2 Qi 2 , (43) 



where 



fin = -2 + ^ {V 2 + ™ 2 - (M 2 - m 2 ) cos 26») - 5e - 2e (e + 2t?" ) , (44) 
fi 22 = -2 + ^2 {v 2 -m 2 + {y 2 +m 2 ) cos 20) +e, (45) 

2 2 

Oi2 = ^21 = - 2H ™ sin 29 - 2e^ (46) 

These equations are solved numerically with initial condition Bn = B 22 ~ k/h, B\ 2 = 0. Figure 1 shows the evolution 
of 



1, 4C 22 1 
sir = — in — — = — m 
2 Cn 2 



^[BMBn] - 5R[S 12 ] 2 ' 



which controls the amount of decoherence, for the two mass ratios /i/m = 1.05 (blue, solid line) and /i/m = 10 (black, 
solid line). This ratio is also the entropy per mode, see section V. Horizon exit occurs at 5 and 3 efolds respectively 



3 For the behavior of real and imaginary parts relating to inflationary perturbations see Q and the next section 
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FIG. 1: The entanglement entropy per mode Sk plotted for the cases fj, = 1.05m (numeric: solid, blue line; analytic: dashed- 
dotted, red line) and = 10m (numeric: solid, black line). For the first case, the analytic and numerical results are hardly 
distinguishable. Horizon exit happens at 5 and 3 e-folds respectively for the two cases and a few e-folds later the density matrix 
has decohered. For similar masses the entropy grows linearly with the number of e-folds until the end of inflation. For dissimilar 
masses, the entropy stops growing when the heavy field starts oscillating and remains constant on average, becoming zero for 
negligible time intervals. 

for the chosen mode and inflation finishes at 60 efolds. We see that Sk becomes larger than unity a few efolds after 
horizon exit and continues to grow linearly with the number of e-fold causing the density matrix to decohere. For 
fj,/m = 1.05 this growth continues until the end of inflation. For the larger mass ratio it stops once the heavier field 
drops out of slow roll and starts oscillating, but s retains its average value until the end of inflation (it becomes zero 
for very short intervals and is quickly restored). Thus, in both cases decoherence is established a few efolds after 
horizon crossing and is retained to the end of inflation. 

The period of oscillations in figure 1 represents the Poincare recurrence time for our system. The recurrence time 
is short because for each mode we have only two degrees of freedom. In models with several scalar fields however, we 
expect that the Poincare recurrence time grows exponentially. Hence the brief recoveries of coherence we see in figure 
1 should be a specific feature of the two scalar field model under consideration. 

The case of almost equal masses can also be treated analytically where both fields exhibit slow-roll: 



d(j) 
dN 



d v V 
3H 2 



dN 



d x V 
3H 2 



1 



3A/2 



-V 



(48) 



with dN = —Hdt = —aHdri, M 2 — (8irG) , H = a/a the (almost constant) inflationary expansion rate of the 
universe. Denoting the number of e-folds N as 



N t -N 



a(N) = a+e 
we have that 

<p = 2M p \/Ncosuj 
X = 2M p VNsmuj. 

As N decreases inflation proceeds, ending when N ~ 1. The anglular variable oj evolves according to 

duj /i 2 — m 2 , 



dN m 2 



■ sin 2to . 



(49) 

(50) 
(51) 

(52) 



Slow roll implies that 



tan< 



• tan bo , 



(53) 
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which gives 

dd fi 2 



■sin 26. (54) 



dN 6H 2 

The above expression can be easily integrated (see [23[ for details). The expansion rate H is then given by 

H 2 = -N (/i 2 + m 2 - (/i 2 - m 2 ) cos 2uS) . (55) 

O 

Furthermore, we are going to need an expression for the slow roll parameter 

1 

>L v 

We have 

,4 , ,,4 I, A i 



^ /H2 = 2^ [^/ dN ^ + ( d x/dNf \ . (56) 



1 m + /i 4 — (/i — m ) cos 2lj 
e = ft 

N [n 2 + m 2 — (p 2 — m 2 ) cos 2w] 

Let us first focus on subhorizon scales, where the d n 6 terms are not important; we have 



(57) 



ihd n B lx = h 2 (B 2 ! + B 2 2 ) - k 2 , (58) 



i 

ihdr,B 12 = h 2 2B 12 B n - {aH) 2 Q 12 , (59) 

with solution Bu = B22 — k/h, B\ 2 = 0, which is the standard Minkowski vacuum for two decoupled degrees of 
freedom. 

On superhorizon scales, an analytic solution to (|4"Tj) - ([43]) can be obtained perturbatively if d n 9 is treated as a 
small parameter. This will be true if 

2 2 

Li 171 , , 

5— < 1. 60 

To 0" 1 order in d v 9 we have 

iftS^Bn ~ h 2 (B 2 n + B 2 2 ) - (k 2 + (aff) 2 Q n ) , (61) 

ihd v B 22 ~ ?i 2 (B 2 2 + B 2 2 ) - (A: 2 + (aH) 2 Q 22 ) , (62) 
i^,Bi2 ^ ft 2 B 12 (B n + B 22 ) • (63) 

Thus, Bi 2 = and -Bn, Bi 2 can be obtained from the single field results (see {H])) where B(k, 77) = A(k,?7)/a 2 and 
in the long wavelength limit 

lim qh ~?——-T{y) . (64) 



Using these, we find to leading order for the correlator 



aiJ / 2tt f k \ 2l/1 . /3 



622 = f frw(^)" 2+ H2^ 2 "' (66) 



where 



^ " 2 + 3e-^2-(M 2 +™ 2 -(M 2 -^ 2 )cos20 o ) (67) 
^2 =s \ - \e - (/i 2 - m 2 + ( M 2 + m 2 ) cos 20 o ) , (68) 
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FIG. 2: Left: An illustration of the Wigner function before decoherence for a highly squeezed state describing long wavelength 
cosmological quantum fluctuations. The field amplitude (horizontal axis) is highly correlated with the field momentum (vertical 
axis). The area of the ellipse corresponds to minimum uncertainty and the entropy is zero. Right: The Wigner function after 
decoherence. The ellipse has increased its area leading to a non- vanishing value for the entropy. Note that the correlation 
between field amplitude and momentum has become less sharp but since An(deco)/' !r cia.Bs "C 1 even after decoherence, the 
prediction for acoustic oscillations in the CMB is not spoiled. 



with tan 0o = ~r tanwo. Then, to first order in d n 9 we have 

i7id„Bi2 ~ ft 2 Bi2 (Bn + B 22 ) + [ih(B 2 2 - B u ) +3aH]d v 0. 
Soon after horizon crossing, Bn and B22 become imaginary and B12 is given by 

N 



B, 



e Js<"™ (Z-f3)aH - 



6H 2 



sin 26> dN, 



where we have taken TV* to be a time when the real parts of Bn and B22 have become negligible, and 

4 [i 2 + m 2 cos 26*o 7 
a = 3 £ 3H 2 = iV' 

u 2 cos29 — m 2 „„ 5 
= 3^ cos20 o = -. 

Using (I55p and (|57|) , and neglecting the dependence of 9 on N in (f70|) we find 

sin 29q 



where 



B 12 (iV) 



I(N) A " J ( :-! 

JV 



1 2 

t fi — m 



h 2\/3 /i 2 + m 2 ~ (p 2 — to 2 ) cos 2luq 



a(N)I(N) , 



5\ e 



N-s 



SI S T + 5 



ds 



1 



1 



1 



1 



= /vv v ur(---7 > jv*)-3r(--7 j jv*)-«r(---7,jv) + 3r(-- 7 ,jv) 

As seen in figure 1, the above formulae reproduce the relevant numerical results with high accuracy. 



(69) 



(70) 



(71) 
(72) 

(73) 



(74) 



V. ENTANGLEMENT ENTROPY 



We now discuss how decoherence can be quantified by the amount of entanglement entropy for cosmological per- 
turbations which is generated after tracing out one of the degrees of freedom. The analysis in [24|j [25| implies that the 
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entropy of a gaussian quantum state can be expressed in terms of the product of the field amplitude and momentum 
variances (A^, A 2 ) as 



S = trln[2A q A ff /n] , 



(75) 



with S = for any minimum uncertainty state such as the state of cosmological perturbations. The variance of the 
momentum A T is calculated from the Wigner function W(q, tt) corresponding to the density matrix: 



W(q,w) = / dndix* pe 



— iAir* -iA* 



which for a gaussian state reads 



W(q, tt) cx exp ( - q-r^q* - (tt - tt c i) (tt - n ct )* ) 



(76) 



(77) 



The classical momentum associated with a value for qi is 



Trie/ = — qiCi2?i . 



(78) 



Note that the existence of a classical momentum, corresponding to a Wigner ellipse rotated in phase space (figure 2), 
does not alter the variances of q and tt and equation ([75]) always holds as can be checked explicitly. 

Before decoherence, A^ = l/3f?[Bn] and A^ = ft 2 5R[Bn]/4 such that the state (fT5"]) corresponds to the minimum 
uncertainty state with vanishing entropy. Calculating the entanglement entropy corresponding to the reduced density 
matrix ([37]), we get for the decohered adiabatic perturbation 



S = V 



d 3 k 
(2^)3 



Sk = V 



d 3 k 1 
In 

(2tt) 3 2 



4C 2 



(79) 



Thus, the entropy is characterized by the same quantity (|47p which controls decoherence and figure 1 also shows the 
evolution of the entropy for the cases of the two mass ratios. 

Equation (|T9"|) establishes a link between the entropy of cosmological perturbations and decoherence; the entropy ([79"]) 
corresponds to the entanglement entropy which equals (minus) the information stored in correlations between the 
adiabatic and isocurvature perturbations. Therefore, it conforms with the standard definition of the entropy of 
quantum systems: tracing out the isocurvature perturbation (the unobservable "environment") generates the entropy 
of the adiabatic perturbation (the system) . 

The effect of decoherence on the Wigner function is illustrated in figure 1. As mentioned above, prior to decoherence 
the state is described by the squeezed ellipse, the area of which is dictated by the minimum uncertainty. After 
decoherence, the variance A qi remains to a good accuracy unchanged but the variance in the momentum has increased 
to A^/^-j = 2C22- This is opposite to the standard treatment of inflationary perturbations where the extremely 
squeezed nature of the state is taken to imply an exact classical correlation between the mode amplitude, drawn from 
a gaussian distribution, and the momentum. Decoherence implies that the momentum should also be drawn from a 
distribution with a dispersion larger than that dictated by the minimum quantum uncertainty. 

The existence of acoustic peaks in the CMB requires that the field amplitude and the corresponding momentum 
exhibit a high degree of classical correlation as was already pointed out in 2Q\. Assuming qi ~ Aqi, we find that 
the existence of such a classical correlation would require the variance of the momentum to be much smaller than the 
corresponding classical value, i.e. 



^7ri (dec) 

kid I 



\/CllC2 



IC 



< 1 



12 



(80) 



This condition is satisfied for both cases as can be seen in figure 3. The second case exhibits an almost perfect degree 
of classical correlation by the time inflation ends and the oscillations of the heavier field have decayed. This behavior 
parallels that of the single field case, where the deviation from a perfect classical correlation of the field amplitude 
and the corresponding momentum is exponentially small The first case also exhibits correlation between field 
amplitude and momentum but which is significantly lower; the variance in the momentum is about 10~ 4 of the classical 
momentum value. In principle, such an effect will tend to wash out the CMB acoustic peaks slightly. Although it 
seems unlikely that such a small relative variance in the momentum would have an observable effect on the CMB, it 
may be interesting to examine its consequences. 
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VI. DISCUSSION 



Let us close with a brief summary of the arguments presented in this paper. We discussed the possible relevance 
of isocurvature perturbations for the decoherence of cosmological curvature perturbations generated during inflation. 
We have argued that qi - $ - q2 can play for cosmological perturbations the role of the triptych "system" - "measuring 
apparatus" - "environment" which is central to considerations of decoherence. The crucial ingredient is the appearance 
of a non-zero imaginary component B12 on superhorizon scales which results in C11/C22 <C 1 exponentially. We 
expect this phenomenon to be generic in all multi-field models since, in general, the imaginary parts are exponentially 
dominant over the real parts. We explicitly checked a free two-field model where no explicit couplings between the 
two fields were postulated. Nevertheless, adiabatic and isocurvature perturbations are generically coupled through 
gravitational interactions, implying that decoherence of this sort will be relevant in general. Since this decoherence 
is very effective, we do not expect the inclusion of non-linear interactions to affect this result. Furthermore, even if 
the entropy - and the momentum variance - increases, the strong classical correlation between field amplitude and 
momentum, responsible for the appearance of the CMB acoustic peaks, persists. 
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